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Abstract. New images of the radio Supernova Remnant (SNR) G65. 1+0.6 are presented, based on the 408 MHz 
and 1420 MHz continuum emission and the Hl-line emission data of the Canadian Galactic Plane Survey (CGPS). 
A large shell-like structure seen in the 2695 MHz Effelsberg map appears to have nonthermal spectral index. HI 
observations show structures associated with the SNR G65. 1+0.6 in the radial velocity range of -20 to -26 km/s 
and suggest a distance of 9.2 kpc for the SNR. The estimated Sedov age for G65.1+0.6 is 4 - 14xl0*yr. The 
pulsar (PSR) J1957+2831 is possibly associated with G65.1+0.6, with consistent distance and kinematic age 
estimate, but different characteristic age than the SNR. The EGRET source 3EG J 1958+2909 and 7-ray source 
2CG 065+00 are also near the eastern edge of the SNR but do not agree in position with the pulsar and are likely 
not associated with the SNR. The SNR's flux densities at 408 MHz (8.6+0.8 Jy), 1420 MHz (4.9+0.5 Jy) and 2695 
MHz (3.3+0.5 Jy) have been corrected for flux densities from compact sources within the SNR. The integrated 
flux density based spectral index (Si/OC!/~°) between 1420 MHz and 408 MHz is 0.45+0.11 and agrees with the 
T-T plot spectral index of 0.34 +0.20. The nearby SNR DA495 has a T-T plot spectral index of 0.50+0.01. 
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1. Introduction 

Since G65. 1+0.6 was discovered as a low surface bright- 
ness SNR (Landecker et al., 1990), little research on this 
SNR has been done. Due to lack of independent distance 
and age estimates for G65.1+0.6, it was difficult for pre- 
vious authors to compare the SNR with nearby discov- 
ered pulsars (Lorimer et al., 1998). The detections of an 
EGRET source 3EG J1958+2909 (Hartman et al., 1999) 
and a TeV 7-ray source 2CG 065+00 (Alexandreas et al., 
1991; Swanenburg et al., 1981) which partly overlap the 
SNR G65. 1+0.6 give increased reason to further study 
G65. 1+0.6. In this paper, we present the SNR's continuum 
images at higher sensitivity and resolution than previously 
at 408 MHz and 1420 MHz, and first investigate Hl-hne 
emission at various radial velocities for detecting interac- 
tions of the remnant with the surrounding gas, estimate 
its distance and age, and access its possible association 
with the nearby EGRET source, 7-ray source and pulsar. 
We see evidence for a large SNR-like shell in the 2695 MHz 
Effelsberg map, and investigate the evidence below. 
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2. Observations and Analysis 



The continuum and HI emission data sets come from 
the CGPS, which is described in detail by Taylor et al. 
(2003). The data sets are mainly based on observations 
from the Synthesis Telescope (ST) of the Dominion Radio 
Astrophysical Observatory (DRAG). The spatial resolu- 
tion of the continuum images of G65. 1+0.6 is 0.8' x 
1.7' at 1420 MHz and 2.8'x5.9' at 408 MHz. The syn- 
thesized beam for the HI line images is 58" x 2.0 ' and 
the radial velocity resolution is 1.32 km/s. DRAG ST 
observations are not sensitive to structures larger than 
an angular size scale of about 3.3° at 408 MHz and 56' 
at 1420 MHz. Thus the GGPS includes data from the 
408 MHz all-sky survey of Haslam et al (1982), sensi- 
tive to structure greater than 51', and the Effelsberg 
1.4 GHz Galactic plane survey of Reich et al. (1990, 
1997), sensitive to structure with resolution 9.4' for large 
scale emission (the single-dish data are freely available by 
|http://www.mpifr-bonn.mpg.de/survey.htmI||). The low- 
order spacing HI data is from the single-antenna survey of 
the GGPS area (Higgs & Tapping 2000) with resolution 
of 36'. See Taylor et al. (2003) for detail of the method of 
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Table 2. 408-1420 MHz T-T plot spectral indices with 
and without Compact Sources(CS) 



64.60 64.50 64.40 

Galactic Longitude 

Fig. 2. Image of box C at 1420 MHz 

combining the synthesis telescopes and single dish obser- 
vations. 

We analyze the continuum and HI images of G65.1+0.6 
and determine its flux densities and distance. For 
G65. 1+0.6, integrated flux density's errors arc found by 
comparing results for several different choices of back- 
ground region. For compact sources, the flux density's 
errors arc taken as the formal Gaussian flt errors. The 
influence of compact sources within the SNR is much re- 
duced by employing similar methods to Tian and Leahy 
(2005). 

3. Results 

3.1. Continuum Emission 

The CGPS continuum images at 408 MHz and 1420 MHz 

are shown in the upper left and right panels of Fig. 1. The 
6 boxes labeled with letters A-F and used for T-T plots 
are shown in the left. The upper right of Fig. 1 shows posi- 
tions and error circles of EGRET source 3EG J1958+2909 
(marked by number 1) and 7-ray 2CG 065+00 (marked by 
number 2). Pulsar J1957+2831 is marked by letter X. The 
lower left shows the 1420 MHz map convolved to the same 
resolution as the 408 MHz map. The 2695 MHz Effelsberg 
map is reproduced in the lower right for reference (Fiirst 
ct al., 1990). The Effelsberg map has a resolution of 4.3' 
and a sensitivity of 50 mKTs- 

The 408 MHz image for G65. 1+0.6 reveals more struc- 
tures in comparison with Landccker's (1990) 408 MHz 
image, especially, side-lobe effects which are obvious in 
Landecker's maps are much reduced on our image. The 
1420 MHz map shows much better detail of the fine struc- 
ture in G65. 1+0.6 than any previous image. More compact 
sources within the SNR are resolved in the 1420 MHz map. 
More filament structures appear in the northern, north- 
western and southwestern parts of the SNR. Fig. 2 shows 
that a ring-shape structure inside box C is clearly detected 



Sp. Index 



Area 


including CS 


CS removed 


A 


0.30+0.21 


0.25+0.22 


B 


0.43+0.21 


0.34+0.20 


C 


0.35+0.13 


0.35+0.13 


D 


0.59+0.13 


0.51+0.20 


E 


0.34+0.59 


0.18+1.20* 


F 


0.50+0.01 


0.50+0.01 


Manual fitt 


ing gives a more reliable value 0.7. 



Table 3. Integrated flux densities and 408-1420 MHz 
spectral indices of G65. 1+0.6 and compact sources within 
the SNR 



Freq. 


G65.1+0.6 


CS of G65.1+0.6 


G65. 1+0.6 and CS 


MHz 


Jy 


Jy 


Jy 


408 


8.6+0.8 


4.3+0.2 


12.9±1.0 


1420 


4.9+0.5 


1.6+0.1 


6.5+0.6 


2695 


3.3+0.5 


0.9+0.1 


4.2+0.4 


a 


0.45+0.11 


0.79+0.06 


0.55+0.10 



near the western diffuse emission regions of the SNR. The 
outhnes of G65. 1+0.6 at 408 MHz and 1420 MHz are very 
similar. 

The 3°x3° 2695 MHz Effelsberg map appears to 
show a circular-shell (Radius about l'^) SNR- like emis- 
sion region including the two known SNRs G65. 1+0.6 and 
G65.7+1.2 (DA495). But the CGPS maps at both 408 
MHz and 1420 MHz don't show this circular shell. We 
mark the large emission region (A) and its two inner ar- 
eas (D and E) related with G65. 1+0.6 in order to study 
the emission further. 

3.2. T-T Plot Spectral Indices 

Bright compact sources affect the measured integrated 
flux densities for G65. 1+0.6 and its measured spectral in- 
dex. Thus we correct for the effects of compact sources. 
Table 1 lists properties of the 16 brightest compact sources 
which are detected in box A, and the first 11 of which are 
within G65.1+0.6 at both 408 MHz and 1420 MHz. 

First we discuss spectral indices between 408 MHz and 
1420 MHz based on the T-T plot method (Turtle et al., 
1962). The principle of the T-T plot method is that spec- 
tral indices {T^=Toi^~^) are calculated from a fit of a 
linear relation to the T1-T2 values of all pixels within a 
given map region. Ti is the brightness temperature of a 
map pixel at one frequency and T2 is for the second fre- 
quency. The higher resolution image has been smoothed 
to the lower resolution for the T-T plot comparison. The 
brightness temperature spectral index /3 is derived from 
the slope of the line. The error in spectral index is derived 
from the uncertainly in slope of the line. The flux den- 
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66.5 66.0 65.5 65.0 64.5 f^^.s 66.0 65.5 65.0 

Galactic Longitude 

Galactic Longitude 



Fig. 1. The first row of images shows the CGPS maps at 408 MHz (left) and 1420 MHz (right). The left of the second 
row shows the 1420 MHz map convolved to the same resolution as the 408 MHz map. The right of the second row 
is the 2695 MHz Effelsberg map. The 6 boxes labeled with letters A-F and used for T-T plots are shown in the left. 
The right of the first row shows positions and error circles of EGRET source 3EG J1958+2909 (95% confidence level, 
marked by mmibcr 1) and 7-ray 2CG 065+00 (90% confidence level, marked by number 2). Pulsar J1957+2831 is 
marked by letter X. The direction of North (N) and South (S) is marked on the lower right image. 



sity spectral index a (S^oci'"") is related to /3 by /3=a+2. 
Spectral index refers to flux density spectral index a in 
this paper unless specifically noted otherwise. 

For the T-T plot spectral index analysis, we select a 
single large region (A) including the two SNRs G65.1+0.6 
and DA495, region B including G65. 1+0.6, region F in- 
cluding DA495 alone and other 3 areas (C, D, E) around 
G65. 1+0.6, as shown in Fig. 1. Region A yields the T-T 
plots shown in the upper half of Fig. 3. Region B yields 
the T-T plots shown in the lower half of Fig. 3. 

In the T-T plot analysis, two cases are considered: us- 
ing all pixels including compact sources; and excluding 
compact sources listed in Table 1 from the images. The 
compact sources are bright compared to the SNR emis- 



sion. Since the compact sources have a steeper spectrum 
than the SNR, they are seen in the T-T plot (Fig. 3 left 
panel) as the steeper lines of points extending to higher 
Tb ■ We completely remove regions of pixels including the 
compact sources from the analysis. Each region is taken 
to be a few beamwidths across, so that any contribution 
from the compact source is below 1% of the diffuse SNR 
emission. 

Table 2 lists the results for two cases of analysis: in- 
cluding compact sources and removing compact sources. 
Visual inspection of the T-T plots confirms that the sec- 
ond method produces the most reliable results. The com- 
pact sources' influence on the spectral index calculation is 
obvious in the T-T plots of Fig. 3, and also seen in Table 
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Table 1. List of the 16 brightest compact sources ( the first 11 sources are within G65. 1+0.6) and their integrated 
flux densities within box A 



No. RA(2000) Dcc(2000) GLON GLAT S408mh^ Si420Mif;. Sp. Index 



[h m s] 



deg deg 



mjy 



mjy 



1 
2 
3 
4 
5 
6 
7 
8 
9 
10 



19 
19 
19 
19 
19 
19 
19 
19 
19 
19 



11 19 



53 46.53 

55 3.77 

53 3.10 
52 53.01 
57 2.74 
52 3.14 

56 11.84 
51 57.62 

54 31.32 

57 7.54 
56 51.73 



27 52 56.6 

28 35 43.2 

29 10 30.2 
28 17 43.0 
28 39 52.6 
28 29 46.9 

27 59 46.5 

28 41 48.9 
28 55 59.3 
28 58 25.4 
28 46 24.5 



64.566 
65.323 
65.593 
64.819 
65.608 
64.898 
64.940 
65.060 
65.551 
65.881 
65.680 



0.109 
0.234 
0.910 
0.490 
-0.102 
0.750 
-0.290 
0.870 
0.510 
0.044 
-0.011 



1388 ±42 
1044 ±34 
610 ±19 
267 ±20 
229 ±16 
191 ±37 
147 ± 9 
137 ±30 
127 ±10 
94 ±10 
83 ±13 



533± 17 
318± 10 
261± 9 
111± 4 
37± 3 
31± 7 
61± 6 
44± 3 
49± 3 
46± 2 
60± 4 



0.77 
0.95 
0.68 
0.71 
1.45 
1.46 
0.70 
0.91 
0.76 
0.58 
0.26 



,73 to 0.80) 
,92 to 0.99) 
,65 to 0.72) 
,64 to 0.78) 
,36 to 1.55) 
,12 to 1.80) 
,61 to 0.81) 
,66 to 1.13) 
,69 to 0.84) 
,49 to 0.67) 
,13 to 0.41) 



12 19 

13 19 

14 19 

15 19 

16 19 



56 46.80 

49 27.70 
49 39.44 
55 3.55 
52 22.96 



28 21 33.9 
28 55 57.2 

28 40 27.0 

29 36 6.0 
29 56 1.9 



65.317 

64.984 
64.783 
66.184 
66.170 



-0.211 

1.461 
1.293 
0.755 
1.424 



765 ±25 

290 ±12 
281 ±12 
117 ± 6 
48 ± 6 



577±340 

135± 7 
100± 4 
56± 2 
42± 1 



0.23( 

0.62 

0.83 

0.75 

0.11 



,15 to 0.94) 

,57 to 0.67) 
,78 to 0.87) 
,70 to 0.80) 
,01 to 0.21) 
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□ flux density of CS 
flux density w/o CS 
best-fit power-law 


□ 


□ 



100 1000 



Frequency (MHz) 

Fig. 4. Radio spectrum of G65.1±0.6 with and without 
Compact Sources (CS). The power law model has a best- 
flt spectral index of 0.49 with la uncertainty range of 0.40 
to 0.57. 

2 for areas A, B, D and E. From now on we discuss spec- 
tral indices derived with compact sources removed, unless 
specifled otherwise. 

Region A includes several emission regions. The spec- 
tral index for a box F including DA495 alone is given in 
Table 2 (0.50±0.01). Prom the right panels of Fig. 3, it is 
obvious that the T-T plot of G65. 1+0.6 is only left lower 
part of region A's T-T plot. The upper right part of region 
A's plot is from the SNR DA495. 



3.3. Integrated Flux Densities and Spectral Indices 

We have derived integrated flux densities of G65.1+0.6 
from the 408 MHz and 1420 MHz maps. Values given 
have diffuse background subtracted. The resulting 408 
MHz to 1420 MHz spectral index, using flux densities 
without compact sources, is 0.45±0.11. Table 3 lists the 
flux densities and spectral indices of G65. 1+0.6 and the 
compact sources within G65. 1+0.6. Compact sources con- 
tribute about 33% at both 408 MHz and 25% at 1420 MHz 
to the SNR's flux densities, and have a signiflcant effect 
on the spectral index. It is noted that the SNR spectral 
index derived from integrated flux densities is consistent 
with the SNR spectral index (0.34±0.20) derived by the 
T-T plot method. 

We have calculated total compact source flux densities 

for 2695 MHz frequencies, using the 408-1420 MHz spec- 
tral index upper and lower limits and flux densities from 
Table 1. The compact sources contribute 21% at 2695 MHz 
to the SNR's flux densities. We have recalculated the flux 
density values of G65. 1+0.6 at 2695 MHz by subtract- 
ing the compact source flux density. We note that Fiirst 
et al (1984) gave a flux density based on the 2965 MHz 
Effelsberg image with resolution 4.4 x 4.4 arcmin, but we 
obtain a new value from the Effelsberg 2695 MHz image 
with a little higher resolution 4.3 x 4.3 arcmin: 3.3±0.5 
Jy and use this instead. We fit the resulting fiux density 
values with a power-law to obtain spectral index. Figure 4 
shows the corrected flux densities and the best-fit power- 
law. The best fit spectral index is 0.49 with la uncertainty 
range of 0.40 to 0.57. The fit to G65. 1+0.6 spectrum in- 
cluding compact sources flux density gives a=0.58. The 
steeper a is expected since compact source has the steeper 
spectrum than SNR (see Table 3). 
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TB(14S0) 



TB(1420) 




TB(1420) 




Fig. 3. 408-1420 MHz T-T plots for region A (the upper half) and region B (the lower half). The left plots for maps 
including compact sources; The right plots for compact sources removed from analysis 



3.4. HI Emission 



nels. Each map has superimposed on it contours of con- 
tinuum emission at 1420 MHz chosen to show G65.1+0.6. 



We have searched the CGPS radial velocity range for fea- 
tures in the HI which might relate to the morphology of 
G65. 1-1-0.6. There is emission which is coincident with the 
boundary of G65. 1+0.6 in the velocity range -20 to -26 
km/s, and only in this velocity range. The spatial rela- 
tion with the edge of G65. 1+0.6 indicates that the HI is 
associated with the SNR. These HI observations show a 
good association of HI features with G65. 1+0.6, similar 
to other accepted HI associations with SNRs and with 
similar HI velocity ranges (~10km/s), e.g. HI associated 
with G126.2+1.6 (Tian & Leahy, 2006), with G127.1+0.6 
(Leahy & Tian, 2006) and with DA530 (Landecker et al., 
1999). Fig. 5 shows maps of HI emission in twelve chan- 



3.5. EGRET Source, j-Ray and Pulsar nearby the SNR 

The upper right of Fig. 1 shows positions of the EGRET 

source 3EG J1958+2909 detected by Hartman et al.( 
1999), the TeV 7-ray source 2CG 065+00 detected by 
Alexandreas et al.(1991) and the PSR J1957+2831 discov- 
ered by Lorimcr (1998). The EGRET and 7-ray sources 
overlap each other and the SNR G65. 1+0.6. The EGRET 
and TeV sources are likely the same source with position 
inconsistent with the pulsar's position and arc likely not 
associated with the SNR. The pulsar is nearby at the east- 
ern boundary of the SNR G65.1+0.6. 



6 



Tian & Leahy: The SNR G65. 1+0.6 and its associated pulsar J1957+2831 



Fig. 5. HI emission in the field centered on G65. 1+0.6 from -20 to -26 km/s and -41 to -46 km/s. The radial velocity 
of each map is indicated at its top left corner. The outline of G65. 1+0.6 in 1420 MHz continuum emission is indicated 
by the contour at 7.2 KTs. 



4. Discussion 

4.1. The Distance and Age to G65.1+0.6 

Since the velocity is negative (-20 to -26 km/s) this implies 
it is outside the solar circle and at 1=65.1'', is very far away. 
We use circular galactic rotation velocity Vi{=Vo=220 
km/s. Then Ro/R=0.89 for -22 km/s, 0.87 for -26km/s and 
0.90 for -20km/s. Then we use Ro=8.5kpc, so d=9.2kpc 
(-22 km/s), 9.0kpc (-20km/s) and 9.6kpc (-26km/s). Non- 
circular motions in the HI feature increase the distance 
uncertainty, eg ±5 km/s, increases the distance range to 
8.7 - 10.1 kpc. This means G65. 1+0.6 is very large (an- 
gular diameters 90' by 51', average diameter 70'): yields 
radius 93.5pc (d=9.2kpc). 

The SNR is of low surface brightness and large ra- 
dius, so probably exploded in a low density part of the 
ISM. A Sedov model (Cox D., 1972) gives t=14xlOV 
if no=0.01cm-3, or t=4.4xlOV if no=0.001cm-3. The 
radius at which cooling affects the dynamics of the SNR 



is 68pc if no=0.01cm^'^, or 127pc if no=0.001cm^''. So 
it is possible that the SNR is still in the Sedov phase or 
has already entered the cooling phase depending on the 
density no. The shock temperature is 8xlO'''K if no=0.01 
and 8x lO^K if no=0.001. Since the density is low it would 
be faint in X-rays, and since it is quite distant it would 
be absorbed in X-rays by the intervening ISM. The non- 
detection in the Rosat All-sky Survey of G65. 1+0.6 is con- 
sistent with a large distance to G65. 1+0.6. 

4.2. Possible Association of the SNR G65. 1+0.6 and 
the Pulsar J1957+2831 

Lorimer et al. (1998) discovered PSR J1957+2831 which 
is located at the eastern edge of G65. 1+0.6, and com- 
mented that the pulsar is not likely to be associated with 
G65.1+0.6 due to the large dispersion measure (DM) dis- 
tance (7.0±2.3 kpc) and old characteristic age (1.6 Myr) 
of the pulsar, and their statistical grounds of pulsar-SNR 
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pairs in the case of absence of G65.1+0.6's distance and 
age parameters. But Lorimer et al. (1998) also admitted 
a possibility that the PSR J1957+2831 may indeed be 
younger than its characteristic age because the pulsar has 
a spectral index typical for other young pulsars. If we as- 
sume that the pulsar birth-place is at geometrical centre 
of G65. 1+0.6, and moves away toward east of the SNR at 
typical velocity of 500 km/s (Lyne et al., 1994), the kine- 
matic age of the pulsar approximately is ISxlO^yr and 
consistent with the estimated SNR's above age. Although 
previous studies (Tian & Leahy, 2004; Mereghetti et al., 
2002; Kaspi et al., 2001) have suggested that pulsar char- 
acteristic ages can be poor age estimators for young pul- 
sars, the pulsar J1957+2831 is not young. The pulsar's 
DM distance (7.0±2.3 kpc) and kinematic age agree with 
G65. 1+0.6, and its characteristic disagrees with the age 
of the SNR, so it is possible that the pulsar is associated 
with the SNR. 

4.3. New Supernova Remnant G64.5+0.9? 

The 2695 MHz Effelsberg map appears to show a circular- 
shell SNR-like emission region most of which is included 
in region A shown in Fig. 1. The 408-1420 MHz T-T plot 
spectral index of areas C, D, E in Table 3 shows they are 
probably all nonthermal emission regions (each is about 
2(7 different from thermal spectral index). This supports 
the possibility that there is a large diffuse SNR including 
the two known SNRs. But the HI images don't show any 
association of HI features with the large circular shell. 
Thus current evidence is inconclusive. 

The ring-like structure in box C detected in the 1420 
MHz image (Fig. 2) may be part of the western edge of 
G65. 1+0.6. However, the ring has a good association with 
HI features in the velocity range -41 to -46 km/s as shown 
in Fig. 5. The HI velocity corresponds to a distance of 
llkpc (-43 kms). Now we name the emission structure 
as a possible new SNR G64.5+0.9 which certainly needs 
to be confirmed by future high resolution and sensitivity 
observations. 
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5. Conclusion 

Using the 408 MHz and 1420 MHz continuum emission 
and the Hl-line emission data of the CGPS, physical pa- 
rameters of G65.1+0.6: flux densities, spectral index, dis- 
tance and age are cither corrected or first obtained. Based 
on distance and age of both the PSR J1957+2831 and 
the SNR G65. 1+0.6, we conclude that the pulsar and the 
SNR are possibly associated. Our analysis of G65. 1+0.6 
shows that the integrated flux density based spectral in- 
dex between 1420 MHz and 408 MHz (0.45±0.11), T-T 
plot spectral index (0.34 ±0.20) and the fit spectral in- 
dex (0.49) are consistent. The nearby SNR DA495 has 
a T-T plot spectral index of 0.50±0.01. A possible SNR 
G64.5+0.9 is suggested for its likely nonthermal emission 
feature, shell-type structure and associating with HI emis- 
sion features. 



